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^ ; ABSTRACT 

t/^ , We investigate the age constraints that can be placed on the double pulsar system using 

models for the spin-down of the first-born 22.7-ms pulsar A and the 2.77-s pulsar B 
with characteristic ages of 210 and 50 Myr respectively. Standard models assuming 
dipolar spin-down of both pulsars suggest that the time since the formation of B is 
^ ' ~ 50 Myr, i.e. close to B's characteristic age. However, adopting models which account 

0^ , for the impact of A's relativistic wind on B's spin-down we find that the formation of B 

took place either 80 or 180 Myr ago, depending the interaction mechanism. Formation 
I 80 Myr ago, closer to B's characteristic age, would result in the contribution from 

J0737-3039 to the inferred coalescence rates for double neutron star binaries increasing 
lO ' by 40%. The 180 Myr age is closer to A's characteristic age and would be consistent 

' with the most recent estimates of the coalescence rate. The new age constraints do 

not significantly impact recent estimates of the kick velocity, tilt angle between pre 
and post-supernova orbital planes or pre-supernova mass of B's progenitor. 
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1 INTRODUCTION 

In addition to its use as a laboratory for studying gen- 
eral relativit y and plasma physics, the double puls ar system 
J0737-3039 (|Burgav et al.ll2003l : itroe et al.ll2004l ) provides 
new insights into the evolution of massive binary systems. 
In the standard model for bi nary pulsar formation (for a re- 
cent review, see Ivan den Heuvel .2007 ). double neutron star 
systems are formed from massive binary systems where the 
initially more massive (primary) star evolves off the main se- 
quence and undergoes a supernova explosion to form a neu- 
tron star. During the evolution of the initially less massive 
(secondary) star, the first-born neutron star accretes matter 
and gains angular momentum spinning it up to short periods 
and (through poorly understood processes; IShibazaki et al.l 
1 19891 ) reducing its magnetic field. If the secondary is suffi- 
ciently massive to explode as a supernova, and the binary 
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system survives this explosion, the resulting system is a pair 
of neutron stars: a short-period recycled pulsar from the pri- 
mary and a young "normal" pulsar from the secondary. 

The double pulsar system J0737-3039 where a recycled 
22.7-ms pulsar (hereafter A) is observed in a 2.4-hr orbit 
around a 2.77-s pulsar (hereafter B) presents a new oppor- 
tunity to study this model. We use the current spin pa- 
rameters of the two pulsars, together with models for their 
spin-down evolution, to place constraints on the age of the 
system. The motivation for this work is twofold. Firstly, a 
better constraint on the system age would reduce the un- 
certainties in empirical studies of the rate of binary neu- 
tron star inspirals — one of the key sources fo r gravitational 
wave observatories (see, e.g.. lKim et al ] |2006l . and references 
therein). Secondly, because the age can be used to directly 
compute the post-supernova orbital parameters, we may be 
able to better constrain B's progenitor mass and in turn 
improve our un derstanding of the formation of this unique 
binary system (jPiran fc Shavivllioosi : IWiUems et al. I l2006l : 
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Table 1. Summary of the model assumptions used in this paper. 



1 - Model 1 



Model 


Interaction 


"A 




Torque decay 




modeled? 








1 


no 


0.0—5.0 


1.4—3.0 


none 


2 


no 


3.0 


3.0 


B; 10 Myr 


3 


no 


3.0 


3.0 


B: 100 Myr 


4 


yes 


0.0-5.0 


1.0 


none 


5 


yes 


0.0-5.0 


2.0 


none 



IStairs et al.ll2006|l. A pre li minar y version of these results was 
presented Lorimer et al.l (|2005l ). In this paper, we consider 
models which account for the modification of B's spin by A's 
relativistic wind. Following a brief review of the properties 
of J0737-3039 in Section [2l we derive age constraints based 
on various spin-down models for both pulsars in Section [3] 
and discuss their implications in Section |4] 



2 J0737-3039 AND THE STANDARD MODEL 

Applying the binary recycling scenario to the double pulsar 
system, we identify A as the first-born neutron star which 
was spun up (recycled) by the mass accretion process. B 
is then the neutron star formed during the supernova ex- 
plosion of the secondary. Using the observed spin parame- 
ters of both pulsars to estimate their surface magnetic fields 

Bsurf = 3.2 X 10^^ G \/pP, we find Bsurf,A = 6.3 x 10^ G 
and Bsurf.B = 1.2 x 10^^ G. That A's field is some three 
orders of magnitude lower than that of B is in accord with 
the recycling hypothesis — i.e. the weaker field of A is a 
consequence of the recycling process. We note however, that 
due to the interaction of A's wind which penetrates deep 
into B's magnetosphere (Lyne et al. 2004), some care should 
be ta ken when int e rpreti ng the exact value of B's magnetic 
field. |a rons et al] (|2005l ) have proposed a model in which 
A's wind exerts a propellor torque on B which dominates its 
spin-down. In this case, the implied magnetic field strength 
of B is a factor of three lower than the above dipole estimate. 

After the accretion phase, it is assumed that both neu- 
tron stars have been spinning down due to a steady braking 
torque; as such they represent independent clocks measur- 
ing the time since accretion ceased. A straightforward test 
of the prediction is to use the characteristic ages of A and 
B: TA = Pa/(2Pa) and tb = Pb/(2Pb). Lyne et al. (2004) 
find TA = 2.1 X 10* yr and tb = 0.5 x 10* yr. Possible 
explanations for this discrepancy are: (i) the standard evo- 
lutionary scenari o does not apply; (ii) as observed in other 
pulsars (see, e.g. jKramer et al. 2003) characteristic ages are 
not reliable due to their simplifying assumptions; (iii) both 
the model and the characteristic ages are wrong. 

Given the circumstantial evidence in favour of the recy- 
cling hypothesis, and the absence of viable alternative mod- 
els, the simplest resolution is option (ii). In the rest of this 
paper, we investigate the consequences of this case and show 
that, when the simplifying assumptions of the characteristic 
ages are taken into account, the apparent age differences of 
the two pulsars can be reconciled. 
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Figure 1. Probability density functions showing the post- 
accretion spin period of A (left panels) and the time since the 
end of the spin-up phase (right panels). For each model, we as- 
sume a range of initial spin periods for B to be anywhere between 
zero and B's current period (dotted curves), and between zero 
and 150 ms (solid curves). 



3 MODELING THE SPIN EVOLUTION AND 
AGES OF A AND B 

Our goal is to use the observed parameters and models for 
the spin-down of the two neutron stars in J0737-3039 to 
place constraints on the system age and pulsar birth param- 
eters. Before describing specifics, we first outline our general 
approach. For each pulsar, we consider a generic spin-down 
model of the form 



(1) 



where P is the spin period, n is the braking index (for spin- 
down due to magnetic dipole radiation, n = 3) and the factor 
K depends on the neutron star moment of inertia, braking 
torque applied to the star and, for some models of B's spin 
evolution, the effect of A's relativistic wind. To distinguish 
between each pulsar, we add A and B subscripts where ap- 
propriate. Our basic approach is to apply and solve the equa- 
tion under certain model constraints to find the 'spin-down' 
age of each pulsar, tsd.A and tsd.B- Table [1] summarizes the 
different models we investigated. 

A key assertion we then make is that the spin-down age 
for A, tsd.A, refers to the time since spin-up ceased, and is 
essentially the same epoch as B began life as a pulsar. We 
therefore assume 



(2) 



Since this condition is only met for certain sets of birth 
parameters, we can use it to find the most probable age of 
the system assuming a given set of model assumptions. Our 
results are summarized in Fig. [T] and discussed in detail in 
the following subsections. 
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3.1 Constant braking parameters 

In the simplest case, wiiere K and n are independent of time, 
equation ^ can be integrated directly to find tsd- For the 
case n = 1, we find 



tsd = 2r In 



while, for all other values of n, the solution is 



2r 



(n-1) 



(f) 



(3) 



(4) 



Here, Po is the initial spin period and the characteristic age 
T = P/2P. Since the current period and period derivative 
are readily observable through timing observations, the un- 
known parameters are n and Po for each pulsar. 

We explore the parameter space using a Monte Carlo 
simulation to compute the spin-down age for B assuming 
an initial spin period and braking index. Then, assuming a 
braking index for A, we compute the required initial spin pe- 
riod of A and the age of the system by asserting equation ([2]) . 
As a starting point, hereafter known as model 1, we adopt 
a braking index for B based on observations of other nor- 
mal pulsars (see, e.g.. iKaspi fc H clfand 2002, for a review) 
which are consistent with a flat distribution in the range 
1.4 < nB < 3.0. Given the completely unknown braking in- 
dices for recycled pulsars in general, for A we took a more 
conservative approach and assumed a flat distribution in the 
range < ua < 5. To show the effect of B's unknown ini- 
tial spin period on the results, we performed all simulations 
using a flat distribution in the range < Po,b < Pb (dotted 
lines) and for < Pq.b < 150 ms (solid lines). The upper 
bound in the latter case is taken from the range of initial spin 
periods i nferred from pulsars with multiple age co nstraints 
(see, e.g.. lMigliazzo et al.ll2002l : lKramer et al.ll2003l '). The re- 
sulting initial spin period distribution for A peaks just below 
20 ms and the age distribution peaks at ~ 50 Myr (i.e. B's 
characteristic age). 



3.2 Exponential torque decay 

Model 1 assumes no decay of the braking torque. An al- 
ternative solution to equation ([1} can be found for the 
case where the braking torque K decays with time. As- 
suming, for simplicity, an exponential decay of the form 
K = A'oexp(-t/idccay), whcrc Ko is a constant and idccay 
is the 1/e decay time, equation ([T]) integrates to yield the 
so-called "reduced age", i.e.: 



^reduced ^dccay lufl -|- ^sd /^dc 



(5) 



Here, t^d is described by either equation or Q depending 
on the assumed value of n. 

The cause and even existence of torque decay in iso- 
lated non-recycle d neutron stars is uncert ain and contro- 
versial (see, e.g., iBhattacharva et al.lll992t ). In young neu- 
tron stars, it is thought to be due to either the decay 
of the magnetic fleld and/or the alignm ent of the spin 
and magnetic axes with time (see, e.g., iTauris fc Konarl 
I2OOII ). Since torque decay is not thought to be significant 
for recycled pulsars after the accret ion phase (see, e.g., 
IBhattacharva k, van den Heuvellll99ll ) , we consider only the 
case in which the torque on B decays. In models 2 and 3 



the simulated distributions shown in Fig. 1 result from the 
equality tsd.A ~ trcduccd,B assuming pure magnetic dipole 
braking (jia = tib = 3) and a torque decay on B with a 
timescale of 10 Myr in model 2 and 100 Myr in model 3. In 
both cases, the time since spin-up ceased is smaller than for 
model 1, with the age distribution peaking at 20-30 Myr. 



3.3 Interaction models 

So far we have assumed the spin-do wn of both pulsars to be 
independent. In reality, as noted bv lLvne et al.l (|2004l ). A's 
rate of loss of spin-down energy is 3000 times that of B; this, 
together with the close proximity of the two pulsars in their 
orbit, means that B's spin-down is significantly affected by 
A's relativistic wind. Direct obs ervational evidence for such 
an interaction was presented bv iMcLaughlin et al. 



(I2OO4I) 



T o model these effects, we follow the results of Lvutikovl 
(|2004l ) and consider two cases. In the first, hereafter model 
4, it is assumed that all the Poynting flux from B is dissi- 
pated when it reaches the interface between A's win d and B's 
truncated magnetosphere. Using equation (9) from lLvutikovl 
(l2004h for this case, we flnd that 



1/3 



(6) 



where fci depends on B's radius and intrinsic magnetic field, 
E is the spin-down energy loss rate of A and D is the sep- 
aration of the two pulsars. The basic spin evolution of B in 
this model implies a braking index tib = 1 which is modified 
by A's spin-down energy lo ss. Such a dependence can also 

be found from the model of lArons et al. (2005') . 

In the second case put forward by Lvutikovl l|2004h , here- 
after model 5, the interface is partially resistive and large 
surface currents produced combine with a poloidal magnetic 
field of B to produce a spin-down torque. This process re- 
sults in a relationship of the form 



Pb = fca 



D2 



1/2 



(7) 



where k2 also depends on B's radius and magnetic fleld 
strength. In this case, the spin evolution of B is independent 
of its own period (corresponding to a braking index tib = 2) 
and is completely dominated by A's spin-down energy loss. 

To derive age constraints for models 4 and 5, we adopt 
a slightly different approach since equations |6] and [7] can- 
not be integrated analytically. Instead, we solve for the spin 
period evolution of B numerically assuming the parameters 
fci and k2 to be constant. Starting with the currently ob- 
served spin parameters, we step back in time and calculate 
the variation of Ea ~ in^ IaPaPa^ using equation (1) to 
compute Pa and Pa at each epoch. We assume the canonical 
neutron star moment of inertia Ia ~ 10"' ^ kg m^ . Simultane- 
ously, we e valuate D usin g the results of lPeters fc Mathews! 
(WM) and iPeterd (| 19641 '). During the calculation, we also 
keep track of whether A's wind continues to pen etrate B's 
magnetosphere using equations (10) and (12) from lLvutikovl 
((200^j3- At the point when this condition is no longer met. 



^ Note that there is a missing c in the numerator of equation (10) 
of lLvutikovl | |2004) to calculate the magnetic field strength of B. 
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Figure 2. The spin history of B assuming two cases (l eft: model 
4; righ t: model 5) of magnetospheric interaction from iLvutikovl 
B's current spin period and epoch is marked by the solid 
point at the top right. The vertical dotted lines in each case show 
the point at which A's wind no longer penetrates into B's mag- 
notosphere (125 Myr for model 4 and 58 Myr for model 5). 



we follow B's spin evolution using the standard spin-down 
formula given in equation ((T} . Example evolution curves are 
shown in Fig. [2] for the case of underlying dipolar spin-down 
(nA = riB = 3). 

The results of our Monte Carlo simulations for models 
4 and 5 are shown alongside the other models in Fig. [1] 
Unlike the other models, the system age constraints depend 
strongly on B's unknown birth spin period. Assuming B's 
period was < 150 ms at birth, the age distributions peak 
sharply at 180 and 80 Myr respectively for models 4 and 5. 
For a broader range of initial spin periods, the system age 
becomes less well constrained. The reason for this can be 
seen in Fig. [2] which shows the large range of system ages 
possible for a given Po,b. 



4 DISCUSSION 

We have considered a range of spin down models to place 
age constraints in the double pulsar system. A striking result 
of this study is the variety of possible system ages. These 
range from ;S 20 Myr (model 2) to almost 200 Myr (model 
4). Models 1-3, which do not account for the effect of A's 
wind on B, favour significantly smaller ages than models 
4 and 5 which do account for the interaction. All models 
we considered favoured an initial spin period for A that 
is close to its currently observed value, with the peak of 
the distribution in the range 15-20 ms. This range is con- 
sistent with A's initial s pin per iod predicted by accretion 
spin-up models (see, e.g.. lArzoum anian ct al. 1999). Given 
the current evidence in favour of initi al spin periods for nor- 
mal pulsars to be < 150 ms (see, e.g.. iMigliazzo et"al]l2002l : 



iKramer et al. 2003 ), and for interaction between A and B 



( McLaughlin et al. 



2004h , we prefer the constraints provided 



by the solid lines for models 4 and 5 shown in Fig. [T] 

Since the age of J0737-3039 determines its contribution 
to the coalescence rate of neutron star binaries, we briefly re- 
vis it the resul t s of th e calculations most recently carried out 
bv lKim et all (|2006l ) where an age of 230 Myr was assumed. 



Taking into account the additional time to coalescence of 
85 Myr, estimates of the total lifetime of J0737-3039 in this 
paper range betw een 105 Myr and 265 Myr. From Table 1 of 
iKim et alT l|2006l) , we find that the contribution to the global 
merger rate made by J0737-3039 would therefore either in- 
crease by 100% in the youngest case, or drop by 10% in 
the oldest case. In our opinion, the most realistic spin-down 
models we have considered are those which take into account 
the interaction, and assume that B's initial spin period was 
negligible compared to its current value. For model 4, the 
estimated lifetime would be ~ 80 + 85 = 165 Myr, i.e. the 
merger rate contribution would increase by 40% over the 
value found by Kim et al., whereas for model 5, the contri- 
bution does not change significantly. Given the uncertainties 
in merger rate estimates, our results do not change the con- 
clusions that binary neutron star inspirals are unlikely to 
be detectable by the current gravitational wave detectors. 
However, the prospects for detection by future instruments 
such as advanced LIGO are excellent. 

Our results can be used to constrain the post-supernova 
orbital parameters in th e double pulsar syste m. Usi ng the 
formulae given by iPeter s & Mathews^ (| 19631 ') and iPetersI 
(|l964h we find the mean orbital separation and eccentric- 
ity after the formation of B to be respectively 1.0 x 10^ m 
and 0.14 for an age of 180 Myr favoured by model 4. For 
the 80 Myr solution from model 5, the corresponding num- 
bers are 0.9 x 10^ m and 0.11. These constraints in turn re- 
strict the allowed rang es of system parameters at the time of 
the secon d supernova ([Piran fc Shavivll2005l : IWiUems et al.l 
2006; St airs et allbooel T 

We have repeated the analysis of 
Tho rsett et al.l l|2005 ^) who constrain the likely kick veloc- 
ity at the time of the second supernova, Vfc, the tilt an- 
gle betweem the pre and post-supernova orbital planes, 5, 
and the pre-supernova mass of B's progenitor star, m-2.i- 
These simulations were carried out assuming the two differ- 
ent prior distributions of the current (unknown) radia l veloc- 
ity of PSR J0737— 3039, as described in detail bv iStairs et al.l 
(l2006l) . 

We have imposed ou r two possibl e age solutions to the 
simulations described by IStairs et al.i (2006) by restricting 
the age ranges considered to be either 70-90 Myr, or 170- 
190 Myr as oppose d to the range of u p to 100 Myr origi- 
nally considered by IStairs et al.l ll20p6|). In general, we ob- 
tain consistent results to Stairs et al. 1 120061 ) which indicates 
that their work does not critically depend on the system 
age. The only exception is the 70-90 Myr simulation as- 
suming the Gaussian radial velocity distribution, for which 
the age constraints favour slightly lower m2,i, Vk and <5 val- 
ues than the unconstrained case. The new age constraints 
derived here are fully consistent with the idea that B's pro- 
genitor was a low-mass star, and that the system received 
a relatively small impulsive kick at t he time of the second 
supernova (jPodsiadlowski et aLllioOsh . 
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